A ring accelerator for matter-wave solitons. 
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We describe a new technique to accelerate and control a confined cloud of ultracold neutral 
atoms. In our configuration, a Bose-Einstein Condensate trapped in a ring is accelerated by using 
an amplitude-modulated optical potential. The spreading of the atomic cloud is avoided by an 
adequate tuning of the scattering length. We also discuss other variations that offer new possibilities 
of robust control of atomic solitons. Our numerical results show that this cold matter accelerator is 
easily accessible in the frame of current experiments. 
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Introduction.- The control of neutral atomic beams 
is a challenging problem in physics due to its 
potential applications in multiple fields like atom 
interferometry 1], coherent control 2], atom clocks Q or 
quantum informational, between others. The most com- 
mon tool used to exert forces in neutral atoms are radia- 
tive interactions |5] that allow to exchange momentum be- 
tween photons and atoms, as in the case of radiative pres- 
sure. Another good example is the dipole force, which 
appears in a laser intensity gradient and which can be 
interpreted in terms of absorption-stimulated emission 
cycles These processes allow to trap and cool atoms 
and are the basis of fundamental experiments with ul- 
tracold atoms like Bose-Einstein Condensation (BEC) in 
alkalii gases 0. 

In the present work we will describe a simple technique 
that can be used to accelerate an ultracold beam of neu- 
tral atoms stored in a ring reservoir 0. The first of these 
storage devices was the nevatron 0, which consists of a 
simple double wire magnetic structure. After the neva- 
tron, other schemes have been developed. In it was 
proposed a circular magnetic waveguide created by four 
coaxial circular electromagnets. In the condensate 
storage ring was obtained by using the quadrupole field 
created by a four- wire geometry. In all these models, the 
condensate moves in regular cycles with constant angu- 
lar velocity and the cloud spreads along the ring during 
beam propagation. The method we propose overcomes 
the previous problems. It allows to increase the angu- 
lar velocity of a BEC confined in a ring and to avoid the 
spreading of the cloud. Thus, the final result is a compact 
beam of ultracold neutral atoms that can be accelerated 
indefinitely with the only limitation of the lifetime of the 
BEC. 

The trick is two-fold: in first place, we propose the use 
of a static amplitude-modulated optical potential. The 
idea is similar to a previous proposal ^3 showing how to 
accelerate a condensate by employing a moving lattice, 
that transmits momentum to the condensate in packets 
of 2kh (being k the wavevector of the optical lattice). 
In second place, to avoid the spreading of the condensate 
during propagation and to get a suitable control of its mo- 



tion, the scattering length has to be tuned to adequate 
negative values. The final result is a compact atomic 
soliton with increasing velocity in each round. The max- 
imum final speed that can be achieved depends on the 
lifetime of the condensate which limits the number of 
cycles. We have also simulated the effect of frequency- 
modulation (FM) lattices in order to show another phe- 
nomena like controllable localization of the cloud. In 
all the cases, our numerical calculations show that the 
parameters we use fit with usual values of current exper- 
iments. 

Thus, we will first introduce the model and equations 
that we have used in this work, then we will present 
the numerical simulations showing acceleration of a beam 
in one and two-dimensional configurations with AM lat- 
tices. We will also include some analytical calculations 
that corroborate the numerical results. Finally, we will 
show the effect of FM lattices that allow to control the 
position of the condensate. 




FIG. 1: [Color Online] Sketch of the system showing the trap- 
ping ring and the optical lattice (in red) used to accelerate the 
atomic soliton. The dipole trap (in yellow) is used to retain 
the cloud before it is accelerated. 



System and modeling equations.- A sketch of our sys- 
tem is plotted in Fig. ^ It consists of a trapping ring 
illuminated with a periodic optical lattice. The mathe- 
matical description of the dynamics of the cloud is given 
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by a Gross-Pitaevskii (GP) equation of the form 



in— - - — V 2 ^ 
dt 2m 



F(r)^ + /7(f)|^| 2 ^, (1) 



where \I/ is the order parameter, normalized to the num- 
ber N of atoms in the cloud: N = J |^| 2 d 3 r, /7(f) = 
4:7rh 2 a(r)/m characterizes the 2-body interaction deter- 
mined by the value of the scattering length a. In our 
model, we consider a ring potential V(f) of radius ro and 
thickness 2Lj_ which varies azimuthally in the following 
form: 



V(r,0) 



V o f(0) r -L ± <r<r 
elsewhere 



(2a) 



trap of width L (shaded region) and a periodic lattice 
connected to it. To put the cloud into motion from the 
trap of depth Vb it is necessary that the energy of the self- 
interactions overcomes the potential energy of the well. 
This yields to the condition: Na cr ~ — 4L i Vn/lC.[l3|. 
The part of the trapped cloud which overlaps the zone 
with a < feels an effective force due to the nonlinear 
interactions. If the value of a is negative enough, this 
force overcomes the square potential and a robust atomic 
soliton is emitted to the lattice. Thus, a has a double 
effect in the condensate: first it acts as a trigger for the 
motion of the soliton and in second place it avoids the 
spreading of the cloud during propagation 16] allowing 
major times for handling the cloud. 



being 



f(0) = 
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Thus, the trap has an angular dependence with the az- 
imuthal variable as it is sketched in Fig. [I]and explicitly 
plotted in the upper-right picture of Fig. With this 
configuration, the BEC is initially confined by a square 
potential in a region of angular dimension — a < < 
along the ring. The scattering length is negative, al- 
lowing the cloud to form a soliton. The trapping well 
is connected to an optical lattice of variable amplitude 
A(0). 

To put the atomic soliton into motion, the value of 
the scattering length is switched to zero in the confining 
region (indicated in orange color in the upper plots of Fig. 
(2J) and the condensate is emitted to the periodic structure 
(see Fig. ||a, b). The process is the same as in atom 
lasers based on control of the scattering length^. Once 
the condensate has left the trap and circulates through 
the ring, the initial confining potential can be removed 
and a is hold at a negative value along the whole cycle. 

Soliton acceleration.- In Fig. [21 we show a one dimen- 
sional version of our model. In this case, we consider 
that V r ^> Vb and thus the cloud is tightly confined in 
the plane perpendicular to the ring. We have analyzed 
the case of a set of N « 10 5 T Li atoms with a nega- 
tive scattering length. We must stress that our results 
are valid for other atomic species like 85 Rb or 133 Cs pro- 
vided that a is tuned to negative values. In our simula- 
tions we have used a = — 2nm, which is easily accessible 
by adequate tuning of Feshbach resonances [l4j. We have 
chosen a ring such that the transverse width of the cloud 
is w± = 0.6/im. The longitudinal dimension of the cloud 
is w ~ 10w±, thus for this cigar shape configurations, 
transverse effects can be neglected in a first approxima- 
tion. 

As it can be seen in the upper plots of Fig. [21 the 
potential along the ring consists of two parts: a square 




FIG. 2: [Color Online] One-dimensional simulation showing 
the acceleration of an atomic soliton. The cloud is emitted 
from the trap at t = Os when the scattering length is switched 
to zero in the shaded region of the upper plots (between dot- 
ted lines in lower images). In a) the amplitude of the lattice 
is constant, yielding to uniform motion. In b) the lattice is 
modulated with an exponential growth, obtaining the acceler- 
ation of the BEC. The horizontal axis corresponds to a path 
I = 200L±. The width of the cloud is w « L = 20L±. In 
the bottom pictures, the vertical axis is time from t = to 
t=1.5s. 

Once the condensate has been emitted from the trap, 
it is necessary that the amplitude A(0) of the lattice in- 
creases along the ring in order to get acceleration of the 
center of mass of the cloud 17]. This can be appreciated 
in the lower pictures of Fig. [21 In Fig. a), we show 
the effect of a lattice with constant amplitude. As it can 
be seen in the caption, no acceleration is obtained. In 
Fig. [21 b) it can be appreciated the acceleration effect 
obtained with an exponentially modulated AM lattice, 
which is a configuration very suitable for experiments. 
We must stress that similar results can be obtained with 
any rapidly growing function, as we have tested numer- 
ically. This acceleration can be explained qualitatively 
with a variational methodfl^. This well-know technique 
uses a Gaussian trial function to average the lagrangian 
density of Eq. As a result, a second order ordinary 
differential equation is found for the center of mass of the 
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condensate, located at an angle 6o(t): 



ArriT'ow ado 

where dots designate time derivatives, m is the atomic 
mass, w is the longitudinal with of the cloud, tq the ra- 
dius of the ring reservoir and I is the following integral: 



exp 



V o f(0)d0. (4) 



Once the condensate has been emitted, if we assume 
that the frequency of the lattice is much faster than 
the oscillations of the condensate, we can approximate 
f(0) = A(6)sin 2 (v6) w A(0)/2 where the factor 1/2 
comes from averaging over the sine function that defines 
the lattice. In fact, our numerical simulations show that 
there is a critical value of the lattice period below which 
the soliton is only affected by the modulation function 
A{0) and not by the lattice. If we assume an exponential 
growth of the form A(0) = exp (rj9), we can obtain a com- 
pact formula for the acceleration of the atomic soliton: 



O = #exp (rjQo) . 



(5) 



being B = r]V e^ w / 2r ^ 2 /(8mr§). Eq. © can be easily 
integrated for initial conditions #o(0) = #o(0) = yield- 
ing to: 



0o(t) = In 



1 - tair \ — - 




(6) 



Eqs. © and © fit with the numerical simulation within 
an error of about 15%. Taking this factor into account, 
the agreement between the numerical simulation and the 
variational calculation is both qualitative and quantita- 
tive. 

We must finally stress that the lattice should be 
switched off and on before the soliton reaches the bound- 
ary in order to avoid reflections. Repeating the process 
periodically, it is possible to obtain acceleration cycles 
that increase the angular velocity of the beam. For the 
experimental parameters given above we obtained a max- 
imum variation of 30mm/ s per cycle of about Is, with- 
out spreading of the cloud and using a static lattice. This 
suppose a a significant breakthrough comparing with pre- 
vious proposals [l£j. 

In order to check the validity of the one-dimensional 
model, we display in Fig. 03 a series of images from a two- 
dimensional simulation in a rectilinear trap. The lattice 
in Fig. 01 has twenty five lines in the space of the initial 
soliton size. For our simulations, this means a fringe size 
of « 0.5/im, which can be easily obtained experimentally. 




FIG. 3: [Color Online] Two-dimensional simulation showing 
the acceleration of a Bose- Einstein soliton in a rectilinear 2-D 
cavity in the presence of an AM optical lattice. Each frame 
has been taken at regular intervals of time from t = to 
t = 2.5s. The numerical values are the same as in Fig. 



The frames have been taken at the same intervals of time 
from t = to t = 2.5s. The other parameters are the 
same as in Fig. |2|b). As it can be appreciated from the 
pictures, the cloud keeps its shape in time, thus consti- 
tuting a two-dimensional Bose-Einstein Soliton. This is 
achieved as a combination of the guiding ring, the lattice 
and the nonlinear interactions. The periodic boundary 
conditions of the simulation allow to repeat the process 
for several cycles. 




FIG. 4: [Color Online] Numerical simulation of the ring con- 
figuration. Each frame has been taken at regular time inter- 
vals from t = to t = 8000Z/J 1 . As it can be seen in the 
pictures, the cloud does not spread during the cycles. 

In Fig. 2] we show the full ring configuration. The 
numerical values are the same as in the previous figures. 
Each frame has been taken at regular time intervals from 
t = to t = 2550^2 1 . As it can be seen in the pictures, 
the cloud does not spread during the cycles. The combi- 
nation of the guiding ring and the nonlinearity forms a 
rotary soliton 20} which can be accelerated in each cycle. 
This opens new experimental possibilities with ultracold 
atom beams like particle colliders or precision atom in- 
terferometers. 

Soliton motion control- The control of the motion of 
the Bose-Einstein Soliton, can be achieved by using AM 
or FM lattices. Fig. [HI shows the effect of an AM lattice 



4 



with linear decrease of the amplitude. Vertical axis is 
time from t = to t = 4500/^_l for (a), and to £ = 
4000/z/j_ for (b). The rest of the parameters are the same 
as in Fig. [3 As it can be appreciated in the picture, the 
soliton turns back due to the decrease in the amplitude 
of the lattice. This effect could have potential use in the 
design of precise atomic interferometers. As it can be 
appreciated by comparing both images a) and b) in Fig. 
the turning point can be precisely adjusted by simply 
changing the frequency of the lattice. 




FIG. 5: [Color Online] Effect of an AM lattice with linear 
decrease of the amplitude. Vertical axis is time from t — to 
t = 4500 /v± for (a), to t = 4000 /u± for (b). The rest of the 
parameters are the same as in Fig. 




FIG. 6: [Color Online] Effect of a FM lattice. Its shape is 
shown for each case at the top of the pictures. Vertical axis 
is time from t = to t = 4000/z/^; where i/± is the radial 
trapping frequency. The rest of the parameters are the same 
as in Fig. 

On the other hand, if the frequency of the lattice de- 
creases along the ring of propagation, the soliton can be 
stopped at a given position. As it is shown in Fig. El this 
allows a robust control of the atomic beam along the ring. 
We have verified that the solution can be generalized for 
a 2-D numerical simulation and different modulations of 
the lattice. 

Conclusions.- In this work we have described a new 
scheme for obtaining acceleration of ultracold atomic 
beams. Our system is capable of both increasing the 
angular velocity of a BEC stored in a ring reservoir while 
avoiding its spreading. We have also discussed other ef- 
fects that can be obtained by using AM and FM lattices 
which allow a robust control of the motion of the conden- 



sate. Our results ate fully testate in the frame of current 
experiments with ultracold atomic gases. 
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